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Abstract 

Sequences of experimental ground-state energies for both odd and even 
A are mapped onto concave patterns cured from convexities due to pair- 
ing and/or shell effects. The same patterns, completed by a list of exci- 
tation energies, give numerical estimates of the grand potential Q(/3, fi) 
for a mixture of nuclei at low or moderate temperatures T — (i~ x and at 
many chemical potentials fi. The average nucleon number (A)(/3, fi) then 
becomes a continuous variable, allowing extrapolations towards nuclear 
masses closer to drip lines. We study the possible concavity of several 
thermodynamical functions, such as the free energy and the average en- 
ergy, as functions of (A). Concavity, which always occurs for the free 
energy and is usually present for the average energy, allows easy interpo- 
lations and extrapolations providing upper and lower bounds, respectively, 
to binding energies. Such bounds define an error bar for the prediction 
of binding energies. Finally we show how concavity and universality are 
related in the theory of the nuclear density functional. 

PACS: 21.10.Dr, 21.10.-k, 21.60.-n, 24.10.Pa 



1 Introduction 

The observation of a valley of stability and the search for mass formulae belong 
to the oldest subjects studied in nuclear physics pQ and continue to be of great 
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interest today [5] . Given the neutron and proton numbers N and Z as indepen- 
dent variables and the corresponding atomic number, A = N + Z, terms such as 
volume energy oc A, surface tension oc As, Coulomb energy oc Z{Z — l)/^! 1 / 3 , 
symmetry energy oc (N — Z) 2 /A, etc., flourish in the literature, and a great deal 
of attention has been dedicated to the consideration of finer corrections, such 
as, for instance, terms s(N, Z) and p(N, Z) that account for shell and pairing 
effects, respectively, and further correlations. This work is motivated by the 
observation that the dominant terms, such as oc A, oc Z(Z — 1) and oc (N — Z) 2 . 
define a notoriously concave energy surface. 

Upper and lower bounds to nuclear binding energies can be deduced from 
such a concavity, provided that deviations from concavity, possibly induced by 
subdominant terms like oc As, s(N, Z), p(N, Z), etc., can be corrected. For the 
sake of simplicity this paper first considers only sequences of isotopes and, thus, 
takes advantage of concavity with respect to N only; Z is frozen. In Section 2 
we begin with a theory at zero temperature and show how elementary, invertible 
transformations of data can generate truly concave patterns. This is obtained by 
an analysis of the table of second differences between binding energies, then by 
a removal of the pairing energy, and finally by an ad hoc, but minimal, parabolic 
term added to the nuclear energies, if necessary. 

Since concavity is also a property of several thcrmodynamical functions, an 
extension of the zero temperature analysis to finite temperatures is in order. In 
Section 3 we discuss properties of that grand potential, Cl(f3,fj,), which can be 
deduced from the experimental data after their tuning. Other thermodynamical 
functions are also considered, and their concavity is tested. Bounds are found, 
and an error bar for predictions is estimated. 

A generalization to concavity with respect to both N and Z is briefly studied 
in Section 4. An additional motivation for our investigation of concavity prop- 
erties is the need, in density functional theories, of concavity, if the universality 
[3] of a density functional must be obtained. A solution will be shown. 

Finally, a discussion and conclusion are given in Section 5. 

2 Concavity for experimental ground-state en- 
ergies 

Our argument is best illustrated numerically, by using a sequence of isotopic 
ground-state binding energies, — We choose to work with the tin isotopes 
(Sn), because they provide a large number of known isotopes for testing our 
method. In addition, in Section 3, we will extend our approach to finite temper- 
atures using a grand potential, which requires that a sufficiently large number of 
excited states are also known. For instance, a table of 25 ground-state bindings 
for the tin isotopes from 110 Sn to 134 Sn reads, in keV, 

{ 934571, 942744, 953531, 961274, 971574, 979120, 988684, 995627, 1004954, 
1011438, 1020546, 1026716, 1035529, 1041475, 1049963, 1055696, 1063889, 
1069439, 1077345, 1082676, 1090293, 1095540, 1102851, 1105320, 1109239}. 
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Despite a standard linear trend because of a "not too much fluctuating average 
energy per nucleon", this list is not making a smooth pattern, even less of a 
concave one. The full line in Fig. 1, where, for graphical convenience, we 
have added to each Ea a constant term, 115 MeV and a linear term, 7.5 x A 
MeV, shows the amount of irregularity in the pattern. The main source of 
irregularity is, obviously, the pairing effect. If it can be removed, concavity 
emerges. Concavity for a sequence of isotones is not a surprise, because of the 
Z(Z—1) nature of the Coulomb term. For a sequence of isotopes, however, there 
is no such obviously quadratic term available, notwithstanding the empirical 
modelization of a symmetry energy oc (iV — Z) 2 . Our choice of Sn, because of 
its long isotope sequence, illustrates this concavity more dramatically. 
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Figure 1: Sn isotope energies: irregular line joining raw experimental ground- 
state energies Ea + 7500 A + 115000; pairing and parabolic corrections give the 
nonconnected dots. 
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Figure 2: Staggering effect, seen from first differences D — Ea+i — Ea + 7500. 

Most energies stated in this paper are in units of keV. Our data come from 
Ref. [1] . See also [SJ ® [7] . These sources usually quote the binding energy per 
nucleon instead of the total binding energy itself and such values per nucleon 
are given to varying numbers of significant figures, from six to seven for the Sn 
isotopes [1] that we consider. Consequently, we quote and use all our binding 
energies to six or seven significant figures. Our values are generally accurate to 
the order of two to three keV near the center of the sequence, and tens of keV 
at both ends for the Sn isotopes. See, in particular, [4 for error bars. 

The staggering effect is clear from the 24 first differences, Ea+i — .EU + 7500, 
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{-673, -3287, -243, -2800, -46, -2064, 557, -1827, 1016, -1608, 1330, -1313, 1554, 
-988, 1767, -693, 1950, -406, 2169, -117, 2253, 189, 5031, 3581}, see Fig. 2. 

The list of 23 second differences (SDs), SD = E A+l - 2E A + E A -i, 
{-2614, 3044, -2557, 2754, -2018, 2621, -2384, 2843, -2624, 2938, -2643, 2867, 
-2542, 2755, -2460, 2643, -2356, 2575, -2286, 2370, -2064, 4842, -1450}, is in- 
sensitive to the constant and linear terms we used for graphical convenience. It 
gives estimates of the "curvatures" of the pattern. It turns out to be far from 
containing only positive numbers. The nonconnected points shown in Fig. 3 
represent this pattern of SDs. A systematic oscillation, reflecting the staggering 
effect, is found. Alternating signs are obviously due to the gains of binding for 
even Sn nuclei because of pairing. The oscillation between SDs centered at odd 
and even nuclei has, roughly speaking, a constant amplitude. Notice, however, 
the maximum in the list, 4842, due to the shell closure at 132 Sn. 
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Figure 3: Second differences SD = E A +i — 2E A + E A _\. Dots, raw data; solid 
line, result of the pairing correction; dashed line, final result after the parabolic 
correction. 



Add now to each even nucleus energy a fixed number, for example p(N, Z) = 
1250 keV, to tentatively suppress the increase of binding due to pairing. The 
resulting list of SDs is attenuated by an amount equal to ±2p, as illustrated by 
the full line in Fig. 3. The attenuated list reads, 

{-114, 544, -57, 254, 482, 121, 116, 343, -124, 438, -143, 367, -42, 255, 40, 143, 
144, 75, 214, -130, 436, 2342, 1050}. All numbers arc now significantly smaller 
than their partners in the previous list of SDs, except for the smaller, but 
still large maximum at 132 Sn. This maximum is positive, however, and causes 
no difficulty for concavity. The interesting point is rather the most negative 
number in the list, namely —143 keV. All negative curvatures can be converted 
into positive ones if we add to every energy an artificial, parabolic correction, 
P x (A - 122) 2 , with P = 75 keV. Incidentally, the lowest point, A Q = 122, of 
the "added parabola" is arbitrary, because SDs will increase by just a constant, 
namely twice the coefficient P of the A 2 term. After this 2P =150 keV shift, 
the whole sequence of SDs becomes positive, 

{36, 694, 93, 404, 632, 271, 266, 493, 26, 588, 7, 517, 108, 405, 190, 293, 294, 
225, 364, 20, 586, 2492, 1200}, see the dashed line shown in Fig. 3. 
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In short, a "concavity ensuring" manipulation for the isotope energies con- 
sists in replacing each energy Ea by E' A = E A + p X Mod[A + 1, 2] + P x (A — 
Aq) 2 + 115000 + 7500 A. (We recall that the terms which we have added are here 
just for graphical, and later, numerical convenience; they do not influence the 
theory.) With p = 1250, P = 75 and Aq — 122, the list of such tuned energies 
E' A (shown by the nonconnected dots in Fig. 1) reads, 

{17479, 13831, 10219, 7301, 4476, 2055, 266, -1252, -2504, -3263, -3996, -4141, 
-4279, -3900, -3413, -2521, -1439, -64, 1605, 3499, 5757, 8035, 10899, 16255, 
22811}. The choice of the two parameters, p=1250 keV and P=75 keV, is 
empirical: one looks for a pairing correction leading to a modest, if not minimal, 
parabolic correction inducing concavity. Other choices for {p, P} are possible, 
but, obviously, within a small range around 1250 and 75. Furthermore, such 
parameters must be readjusted for different regions of the table of nuclei, but 
it is again obvious that readjustments will be moderate; for instance, the order 
of magnitude for pairing will always be around ~ 1.0 to ~ 1.5 MeV. Analyzing 
short sequences obviously leaves fewer negative SDs to be compensated by the 
artificial, parabolic term, and, hence should induce smaller values of P. For this 
reason, we would expect short sequences to give often better extrapolations. 

Once concavity is obtained, it is straightforward that extrapolations from 
two points on the concave pattern allow predictions of lower bounds to nuclear 
energies and interpolations give upper bounds. Then, from such bounds for 
energies E' , one recovers bounds, of strictly the same quality, for the physical 
energies E. This obtains by subtracting from each E' bound its "tuning term". 

Assume now that 110 Sn were unknown and one had done a brute force ex- 
trapolation, 2-Em - Eu2 = -931957. Compared with Euo = -934571, this 
prediction underbinds 110 Sn by 2614 keV. Consider rather 2E{ n -E{ 12 = 17443, 
to be compared with E( 1Q = 17479. A slight overbinding, by 36 keV, is found. 
Although this small error is likely accidental, it is clear that systematic lower 
bounds will be found. Naturally, once a value for E' A is predicted, one recov- 
ers as good an estimate for Ea after removing the "concavity manipulation 
terms" , which are known explicitly. In the present case, the tuning of the data 
added 952050 to E\\q. The same 952050 must be subtracted from that value, 
2E(u — E[ 12 = 17443, extrapolated from the concave pattern, yielding a fi- 
nal result of —934607, to be compared with Euo = —934571, which obviously 
exhibits the same slight overbinding, 36 keV. Based on the data in [4], the exper- 
imentally measured energy for 110 Sn has an error bar ±14 keV; hence, in this 
case, the experimental error bar and the uncertainty in our theoretical value 
have the same order of magnitude. 

If 132 Sn were unknown, the brute force extrapolation gives, 2Ei 3 i — £130 = 
— 1100787, to be compared with E 132 = —1102851, showing an underbinding 
equal to 2064 keV, clearly failing to reproduce the shell closure effect. Now, 
from the concave pattern, we obtain 2E{ 31 — E{ 30 = 103 1 3, to be compared with 
E{ 32 — 108 99, producing an overbinding by 586 keV. The same overbinding 
is obviously found if one subtracts from 2E{ 31 — E[ 30 = 103 1 3, the tuning 
difference, E{ 32 — E i32 — 11137 50, with the result, —1103437, to be compared 
with -E132 = —1102851. See Fig. 4, which also illustrates an extrapolation using 
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i?i33 and £'134 and an interpolation using E131 and -E133. 
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Figure 4: Left, raw energies: full line, interpolation from 131 Sn and 133 Sn; long 
dashes, extrapolation from 130 Sn and 131 Sn; short dashes, extrapolation from 
134 Sn and 133 Sn. Dots represent the experimental energies E A . Right: same 
inter- and extrapolations, concave data. 



Let us now also compare for an odd rather than even isotope, e.g., 117 Sn, 
two extrapolations and one interpolation from raw data and concave data, re- 
spectively. With our choice of parameters, the tuning difference, E( 17 — Em, is 
994375, with E\n — —995627. Extrapolations from lower masses give, 2En§ — 
Ens = -998248, and 2E( 16 - E{ 15 - 994375 = -995898, hence, overbad- 
ings of 2671 and 271, respectively. Extrapolations from higher masses give, 
2Ens - Eua = -998470, and 2E{ 18 - E{ w - 994375 = -996120, leading to 
overbindings of 2843 and 493, respectively. Interpolations from the nearest nu- 
clei give, {E n& + Ei W )/2 = -996819, and (E{ w + E{ ls )/2 - 994375 = -995494, 
with now overbinding of 1192 and underbinding of 133, respectively. The poor 
results from raw data stress the importance of pairing corrections. The much 
better results from the concave pattern do take this pairing into account, but, 
furthermore, they again make it trivial that extrapolations give overbindings, 
while interpolations give underbindings. 

For the sake of completeness, we now consider a case with a two-step ex- 
trapolation. If the tuned pattern followed a perfect parabolic form, E' A = 
a A 2 + b A + c, the formula, E' A+2 — 3 E' A + 2 E' A _ l —6 a, predicts, for the lower 
bound, an error which is a little larger than desirable. This is because of the 
coefficient, 6, and the fact that the "curvature parameter", a, has here an order 
of magnitude in the few hundred keV range, see those residual second differences 
used for the dashed curve in Fig. 3. This indicates that further stages of the 
theory, with polynoms at least, are worthy of consideration. However, to remain 
within the "linear" scope of the present paper, we give the results obtained for 
135 Sn, -1113208, and 136 Sn, -1115027, when we use the lower bound estimates 
from the formula, 3E A — 2E' A _ V From Ref. [1], the predictions, which are 
based on neighbor nuclei with the same Z and N parities, are —1111320 and 
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— 1115064, respectively, with error bars of ~ ±.5 MeV. We do find a lower bound 
for 135 Sn, with apparently a severe, but mostly expected, error of a couple of 
MeVs. For 136 Sn, we apparently fail to generate a lower bound, but by so small 
a difference that our results stay well inside the error bars. 

To summarize Section 2, when concavity occurs, one concludes that extrapo- 
lations from two points on the concave pattern allow predictions of lower bounds 
to nuclear energies. In the same way, interpolations provide upper bounds. The 
quality of such "concavity bounds" from zero temperature data is found to be 
good; interpolations and extrapolations from the raw, nonconcave pattern, are 
clearly less satisfactory. But there is a more profound reason why a concave 
pattern is necessary. Indeed, several thermodynamical functions, governed by 
theorems proving their concavity, have a notoriously singular limit at zero tem- 
perature: they become nonanalytical and are just piecewise continuous. Their 
limit plots are made of segments; derivatives are discontinuous at turning points. 
Because of the staggering effect, the concave envelope of the raw pattern of Ea's 
would contain only the even isotopes. Concavity is, thus, necessary for a theory 
that must accommodate both odd and even nuclei. Therefore, the thermody- 
namical functions studied in Section 3 preferably use concave energies E' A and 
the corresponding excited state energies E nA . 



3 Concavity with thermodynamical functions 

Consider the particle number operator A and a familiar nuclear Hamiltonian 
H = Ylt=i t'i + Si>j=i v ij> where A, t and v are the mass number, one-body 
kinetic energy and two-body interaction, respectively. Nuclear data tables [HI [H] 
give precise values for a large number of lowest-lying eigenvalues E nA of H, for 
many nuclei. One may, thus, reasonably estimate the grand partition function, 

Z(ji, (3) = Tr exp \J3 (//A - H)] = ^(2j nA + 1) exp \fi (pA - E nA )} , (1) 

nA 

provided that i) the temperature, T — /3 _1 , is low enough to allow a truncation 
of the spectrum to include only those states provided by the tables and ii) the 
chemical potential, p, selects mainly those nuclei in which we are interested. Let 
( ) denote, as usual, a statistical average. The (equilibrium!) density operator 
in Fock space, p = Z^ 1 exp \J3 (pA — H)] , ensures that the following grand 
potential, /?) = ( (H — pA) } — TS, is minimal in the space of many-body 
density matrices with unit trace, since, by definition, (A) = TrpA, (H) = 
TrpH, with the entropy, S = — Tr (p logp). 

We shall rather use (— f2) in the following, to make upcoming proofs of 
concavity slightly easier. This grand potential also reads, 

- Sl(ji, (3) = /T 1 In Z = /T 1 In |^(2j„ A + 1) exp {(3 (pA - E nA )\ j , (2) 
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Simple manipulations then give the relevant statistical averages ( ) of particle 
numbers and energies, together with their derivatives and fluctuations, 



d(-Q)/dfi = (A)=J2 A PA, PA = Z- 1 Yp-inA + 1) exp \p (nA - E nA )] , 

A n 

(3) 

9 2 (-0)/V = d(A)/d f , = [3 ( (A 2 ) - (A) 2 ) , (4) 



and 
then 

or as well 



<9(-0)/<9T = S = ln£-/3(GuA-H)), (5) 
(H) = Z- 1 J2(2j nA + l)E nA exp [P (fiA - E nA )} . (6) 



Furthermore, 

d 2 (-n)/dT 2 =[3* ( (fiA - H) z ) - ( (fj,A - H) } 2 , (7) 



and 



d 2 (~n)/(d^ dT) = -/3 2 [( A (fj,A — H) ) - (A) ( OA - H) )] . (8) 



For our investigations, we will freeze f3 as real and consider functions of a real 
fi. It is then well known that (—fl) is a concave function of fi and T. In turn, 
the double Legendre transform, with respect to both \i and T, 

Hd(-n)/dn + Td(-S})/dT + n = (H), (9) 

shows that (H) is a concave function of both (A) and S, the conjugate variables 
of /x and T, respectively. 

In the following, we do not perform the full, double Legendre transform. 
We rather retain an intermediate representation, with (A) and either T or 
f3. We stay with real variables and functions. We stress that, while A has 
a discrete spectrum, (A) is continuous, a monotonically increasing function 
of /i, smooth provided (3 is finite. The monotonicity results from Eq. (|4]). 
Actually, at low temperatures, strong derivatives signal the onset of discrete 
jumps due to the integer spectrum of A, but we may stay away from this 
"jumpy" regime in the following, at least temporarily. Anyhow, at any fixed, 
finite /3, the smoothness and monotonicity of (A) with respect to /i allows a 
reasonably easy numerical calculation of the inverse function ^i((A}) . Thus, a 
main argument of this Section is that at fixed temperatures we will use (A) as 
a continuous variable and attempt extrapolations towards unknown nuclei. 

For this, given a value of T, we keep track of (A) and (H) as functions of 
/i. Since the functional inversion from (A)(/x) to fi((A)) is reasonably easy, we 
can plot (H) in terms of (A) and attempt an extrapolation for further values of 
(A) . This extrapolation can be considered as a "candidate" for a mass formula, 
at that finite temperature T. 
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According to Eq. ([7]), the average, constrained energy, ((H — /xA)), is a 
monotonically decreasing function of 0. Furthermore, at least for negative chem- 
ical potentials /it, and, more generally, if A has an upper bound, the operator, 
H — /iA, is bounded from below. Therefore, there is a convergence of the process 
consisting in i) extrapolating with respect to fi both ((H — //A)} and (A) for 
fixed values of 0, then ii) eliminating /i to generate the /3-parametrized "mass 
formula" (H) ( (A) , 0) , and finally iii) considering the limit of this mass formula 
when — ► +oo. Alternately, it is equivalent, and maybe more efficient, to first 
eliminate n and then extrapolate the "mass formula" (H) ( (A) , 0) , first with 
respect to (A), then with respect to 0. 

Is there concavity in this intermediate representation? Clearly, a simple 
Legendre transform of (— f2), with respect to fi only, returns a free energy, 
F = (H) — T S, as a concave function of (A) and T. If T is low enough to allow the 
product T S to be neglected, then, at fixed T, one may accept that (H) is an "al- 
most" concave function of (A) . This assumption will be tested by the numerical 
results which follow. Incidentally, a straightforward calculation of A" ((A), 0) = 
d 2 {H)/ (9(A)) 2 yields, A" cx ((A 2 ) - (A) 2 ) (A 2 H) + ((A)(A 2 ) - (A 3 )) (AH) + 
((A)(A 3 ) - (A 2 ) 2 ) (H), with a positive factor, ((A 2 ) - (A) 2 )~ 3 . More simply, 

|!M((A>,/3) cx <(AA) 2 )((AA) 2 AH)-((AA) 3 )(AAAH), (10) 

if one uses the centered operators, AA = A (A) and AH = H (H). From 
Eq. (fT0|) . concavity for (H) is unclear; we shall have to test it numerically. 




Figure 5: Energy (H)' (full curves) and free energy F' (dashed ones) as functions 
of (A)', calculated from concavity tuned energies E' nA . Dots represent tuned, 
ground-state energies, Eq A = E A , already used in Fig. 1 and the rhs part of 
Fig. 4. Upper full and lower dashed curves, T = 150 keV. Lower full and upper 
dashed curves, T — 60 keV. Notice how the full curves turn out to be concave. 

The results, which follow, use [H [9] for the excited-state energies of the 
Sn isotopes. Using the first 10 levels of 110 Sn to 132 Sn, hence a maximum 
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excitation of, typically, ~ 1.2 MeV for odd isotopes and twice as much for 
even ones, we calculate Z, see Eq. (fT]). For those rare cases, where the spin 
j n A is ambiguous, we choose the lowest of the suggested spins. If the spin is 
completely unknown, we set it to be either or 1/2, according to A. These 
tactics minimize the statistical influence of such rare cases. The highest excited 
levels that we use have an excitation energy of order a few MeV; hence, we 
restrict our thermodynamical analysis to temperatures of order ~ 500 keV at 
most. It should be noted that we use the concavity tuned energies, E' nA — 
E nA + 115000 + 7500 A + 1250 Mod[A + 1, 2] + 75 (A - 122) 2 . Consequently, 
all our calculated quantities, (A)', S', (H)', etc., derived from Eqs. ©, ©, 
and ©, for instance, carry a prime superscript. We want to stress that the 
spurious exponential factor created by the constant, 115000, appears in the same 
way in numerators and denominators and, hence, cancels out; it only helps the 
graphics and numerics. Similarly, the linear term, 7500 A, means but a change 
of reference for the chemical potential. Only the pairing and parabolic tuning 
terms make the primed quantities truly distinct from those "raw" , unprimed 
ones, obtained with the untuned energies E nA . 

As an example of our results we choose for Fig. 5 the sequence of six iso- 
topes from 119 Sn to 124 Sn, because of its still severely contrasted list of second 
differences, {588, 7, 517, 108}, after tuning. This produces a serious geometrical 
constraint on the plots of (H)' and F', as they reach their zero temperature 
limit, and better validates our approach. The upper curves in Fig. 5 are the 
plots of the function (H)'((A)') when T = 60 keV (lower full curve) and T = 150 
keV (upper full curve), respectively. The increase of (H)', when T increases is 
transparent. The striking result is the apparent concavity of both curves. It 
is also found that the lower temperature, 60 keV, is low enough to allow (H)' 
to run, in practice, almost through the ground-state energies of both even and 
odd nuclei. As a check of our results, we repeated our calculations using the 
untuned energies, i.e., those which lack concavity. We then found that the low 
temperature limit of the (H) curve went through the even nuclei only, namely, 
the limit curve follows the concave envelope of the experimental pattern. 

For graphical and pedagogical convenience, Fig. 5 shows the plots for six 
nuclei only, but the same observations hold for full plots, with 110 < (A)' < 134. 
As a test, we also calculated (H)' and its low temperature limit when the levels 
of 132 Sn are omitted from the trace sum, Eq. |T]). In that case, the concave 
envelope goes through the dots representing the two odd nuclei, 131 Sn and 133 Sn, 
and ignores the dot representing 132 Sn. Similar verifications of other concave 
envelopes were obtained by removing other nuclei. 

To verify whether concave envelopes and concavity for (H)' result from neg- 
ligible values of the entropy term in the free energy, or, more precisely, negligible 
values of its second derivative, a calculation of TS' is in order. Figure 6 shows 
the difference, TS', between the energy and the free energy, as a function of 
(A)', for T = 60 and 150 keV, respectively. At the higher temperature, 150 
keV, wiggling effects seem to be small enough to allow for the concavity of (H)'. 
At the lower temperature, 60 keV, the wiggling is stronger. As a consequence 
of such a strong wiggling, possibly translating into strong second derivatives, it 
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is not excluded that SDs coming from T S' might prevent (H)' from having the 
proven concavity property of the free energy We must, therefore, numerically 
calculate second derivatives, see Eq. (|TU)) . 




Figure 6: Entropy energy TS'((A)'). Upper (lower) curve, T = 150 (60) keV. 

In Fig. 7 we show, for tuned data, plots of the rhs of Eq. (fT0|) , divided by 
((A 2 )' — (A)' 2 ) . This quantity must remain positive if (H)' shows concavity. 
The full line represents the situation when T = 60 keV, the dashed line cor- 
responds to T = 150 keV. We see that the second derivative, <9 2 (H)'/<9(A)' 2 , 
remains positive almost always. For low temperatures, however, negative values 
may appear. For instance, the full curve, corresponding to T = 60 keV, indi- 
cates small, but definitely negative values around 119 Sn. Numerical tests, which 
are not easy because of a difference effect between the two terms present in Eq. 
(fT0|) . show that the occurrence of such negative, actually moderate, values for 
lower temperatures might be somewhat frequent, while not systematic. Further- 
more, such "negativity accidents" turn out to be worse when we use untuned 
data, maybe because the untuned data lack concavity in the first place. 
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Figure 7: (a) Behavior of H" = <9 2 (H}'/d(A)' 2 ; full (dashed) curve, T = 60 
(150)keV. Notice tiny excursions in negative value territory by the low T curve. 

A likely reason for the negativity accidents with tuned data might be that 
the fluctuation of A is not large enough to justify our use of (A)' as a continuous 
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variable. Since it interpolates between integers, a fluctuation of order ~ 1, or at 
least ~ .5, might be necessary. As shown by the plots in Fig. 8, corresponding to 
T = 60 keV and 150 keV for the lower and upper curves, respectively, a minimum 
temperature of order ^150 keV is needed to avoid too small a fluctuation of 
the particle number. Indeed, for, e.g. ~ 100 keV, we found that the fluctuation 
is smaller than .5 for almost half of the calculated values of (A)'. 

At this stage, the situation can be summarized as follows. On the one hand, 
the tuned pattern of experimental energies shows concavity, but the concavity of 
(H)' as a function of (A)' is not sure, although it seems to occur most of the time. 
On the other hand, we have a theorem proving concavity for the free energy, 
either F = (H) -T S = fi (A) + tt or F' = (H)' -T S' = fi (A)' + 0', as functions 
of T and (A) or (A)', respectively. For instance, elementary derivations show 
that, in that representation where (A) (or (A)') and 8 are the primary variables, 

d¥' _ 2 OF' d 2 F' _ T 

dp ~ * * ' d(A)< ~ M and 5(A)' 2 - <A2>' - (A)'2 ' [n) 

To summarize this discussion, we see that the removal of the entropy term, 
leading from the free energy to just the energy, can sometimes destroy the 
concavity depending upon the temperature, but only weakly, e.g., see Fig. 7. 
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Figure 8: Particle number fluctuation. Upper (lower) curve, T = 150 (60) keV. 

A compromise may be found so that T is low enough to make TS' small 
with respect to F' and high enough to ensure both sufficient values of AA and a 
tolerable oscillating of TS'. Our rationale will be that (H)', even though it might 
deviate from concavity, will stay close enough to the concave F'. Their difference, 
TS', a positive quantity, will define an error bar between a lower bound F' and 
an upper bound (H)' for ground-state tuned energies. Consider again Fig. 5. 
There is no need to stress how, given T, the energy and free energy curves 
make a band defining upper and lower bounds for the experimental energies. 
Moreover, a similar, but narrower band is obtained if T decreases. Figs. 5, 6 
and 8 suggests T ~ 150 keV for an error bar ~ 500 keV between energy and 
free energy, a smooth enough entropy energy and a sufficient fluctuation of the 
particle number. This order of magnitude for T is compatible with the average 
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level splitting, of order ~ 250 keV, that is observed from the first 10 levels of 
110 Sn, ... , 132 Sn used in our calculations. 

The following properties, i) the average energy and the free energy are in- 
creasing and decreasing, respectively, functions of increasing T, ii) the energy 
is larger than the free energy and iii) the entropy term by which they differ 
vanishes when T vanishes, are not big surprises. It can be concluded that, in 
so far as thermodynamical functions can be calculated at low enough temper- 
atures, precise "accuracy bands" may be available. Their, hopefully analytic, 
continuation for higher and/or lower values of A than those available in nuclear 
tables provides a prediction scheme for exotic nuclei. 

An estimate of the entropy term is now useful. Given [i and a large /3, let Aq 
and Eq correspond to that nucleus whose ground-state energy maximizes the 
exponential, exp[/3 (/j,A—Eq A )]. Consider now the first subdominant exponential. 
It might be generated by the first excited state of the same nucleus, or by the 
ground-state of one of its neighbors. Let A\ and E{ be its parameters and define 
A = /i (Ai — A ) — (E{ — Eq). Concavity guarantees that A < 0. Whenever 
is large enough, it is trivial to reduce the grand canonical ensemble to a two 
state ensemble, and the entropy then boils down to s = — e /3A /3A. Hence, the 
product, T s — — e 13 A A, vanishes exponentially fast when /3 — > oo. The rate of 
decrease is governed by that scale defined by A, to be extracted from the tuned 
data. Then one can estimate an order of magnitude for the difference between 
the free energy and the energy. This estimate can be viewed as an error bar for 
the prediction of exotic nuclei via the present "concavity method" . 

We have applied the approach to several other isotopic regions, such as 
Sm, Hg, Pb. In all cases, making the switch to the concave shape improves 
our ability to make extrapolations and interpolations for unknown bindings. 
However, estimating the free energy and the average value of H, to set narrow 
upper and lower bounds on these binding energies, remains more difficult in 
drip line regions, because of the lack of experimental data on excited states in 
such regions. Very simple linear extrapolations still provide, nonetheless, error 
bars with an order of magnitude AE ~ ±400 keV, similar to the best "band" 
seen in Fig. 5. A few examples of the kind of results that we can obtain from 
this thermodynamical error-bar analysis, when using the method described in 
Sec. 2, are these conservative estimates of the bindings of m Hg, 208 Hg and 
209 Hg, namely 1314.33 MeV ±400 keV, 1629.13 MeV ±400 keV and 1632.16 
MeV ±400 keV, respectively. 

4 Two-dimensional analysis; density functional 

We now want to consider binding-energy systematics in two dimensions, i.e., 
N and Z. The second differences studied in this Section are, SD n = En+i.z + 
En-\,z — 2Enz, SD P = En.z+i + E^.z-i — 2Enz, SDb = {En+i.z+i + 
£jv-i,z-i)/2 - E NZ and SD S = {E N+ltZ -i + £;v-i,z+i)/2 - E NZ . If N and 
Z were continuous, and, furthermore, if Enz were a smooth function of N and 
Z, the numbers, SD n , SD p , SDb and SD S , might be interpreted as estimates 
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of second derivatives of the function E(N, Z) in the neutron, proton, constant 
N — Z and constant A directions, respectively. Their connections with neutron, 
proton and p-n pairings and with a symmetry energy are also transparent. 

Setting aside what we consider to be "light" nuclei, namely those of the s-, 
p-, sd- and f 7 / 2 shells, we calculated these four SDs for the rest of the known 
nuclear table. Their maxima and minima are found to be of order 
~ {4840, 6180, 4530, 5900, -3670, -3080, -3160, -2150} keV, respectively. The 
lhs parts of Figs. 9-12 show, in scattered plots with respect to A, the patterns 
of those SD values. The negative SDs, contradicting concavity, are due mainly 
to pairing effects. An extension of the correction done in Sec. 2 is in order. Let 
p be a mass dependent [TU] estimate of pairing. Diminish the binding of doubly 
even nuclei by p, increase the binding of doubly odd ones by p, and leave odd 
nuclei untouched. We set, empirically, p — 1160 — A. There is no doubt that 
better parametrizations of p are possible [10] , but the present linear decrease 
suffices for our demonstration. With this correction, the worst SDs become 
~ {2790, 3970, 2470, 3680, -1500, -930, -980, -120} keV, respectively. As ex- 
pected, this is a large reduction of both positive and negative amplitudes. This 
is shown by the rhs parts of Figs. 9-12, respectively. 
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Figure 9: Scatter plot of convex second differences SD n in the neutron direction 
throughout the table of known nuclei for 27 < A < 251. Left: with bare energies. 
Right: after a pairing correction. 
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Figure 10: Same as Fig. 9, but for SD P in the proton direction. 
While the pairing correction has not completely suppressed negative SDs, 
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note the worst minimum is still of order ~ —1500 keV, concavity has been 
improved significantly and extrapolations are likely made more secure. 

Similar to our one-dimensional analysis in Sec. 2, one can now add to the 
energies a paraboloid term, of the form, u (N — No) 2 + 2v (N — No) (Z — Zq) + 
w (Z — Zo) 2 , in order to get rid of the residual, negative SDs. (We recall that 
SDs do not depend on No and Zo, which are only parameters for graphical 
and/or numerical convenience.) Notice from Fig. 12, and the corresponding 
worse minimum, of order ~ — 120 keV, that the direction where A is constant 
does not deviate much from concavity. The dominant term in the paraboloid 
must, therefore, indicate a direction approximately orthogonal, and also take 
into account the bend of the global stability valley towards neutron excess. 
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Figure 11: Same as Fig. 9, but for SDb in the constant N — Z direction. 




Figure 12: Same as Fig. 9, but for SD S in the direction where A is constant. 

There is, furthermore, no doubt that the two-dimensional free energy, now 
a function of both TV and Z made continuous by a grand canonical ensemble, 
remains a concave function, available for lower bounds. This is a theorem, which 
can be proven in a straightforward manner following the previous derivation for 
one dimension. Simultaneously, it is likely that the mean energy will be again 
"almost concave" , to yield upper bounds, in which case the same error bar that 
was estimated in Sec. 3 remains valid. 

A systematic use of such properties throughout the nuclear table is way 
beyond the scope of the present paper. Rather, we want to discuss the relation 
of concavity with the derivation of a nuclear density functional. Recall that, 
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given an arbitrary many-body density matrix M in Fock space, the density 
functional derives from the definition [llj . 

F\p p ,p n ] = InfM^,p„TrHM. (12) 

Here the symbol, M — > p n ,p p , means that a minimization of the energy is 
performed upon density matrices having the same neutron and proton pro- 
files p n ,p p . A nucleus is identified by the integrals, N — J drp n {r) and Z = 
fdrp p (r), and a later minimization is performed with respect to the profiles 
pn, P P under such identification constraints, 

E Z N = Infp„— JV.ft^Z F[Pn, Pp}- (13) 

The functional T should be universal, in the nuclear physics sense. Namely, it 
should not depend on N and Z. (In atomic and molecular physics, universality of 
the DF has a different meaning; it refers to arbitrary external potentials.) While 
BCS and Hartree-Bogoliubov calculations allow a distinction between even and 
odd particle numbers, universality in density profile space does not allow such a 
separate treatment. Observe, furthermore, that, because of the fact that there 
are many convexities in the pattern of raw energies, there exist many cases 
where, for instance, three nuclei, Ax , A 2 , A3, are such that N 2 = {Nx +N 3 )/2 and 
Z 2 = [Zi + Z 3 ) /2 and {Ex + E 3 ) /2 < E 2 . Then, in a search for E 2 , the mixture 
density matrix, M m ix — { \Ax)(Ax\ + \A 3 ){A 3 \) /2, provides the correct average 
particle numbers N 2 ,Z 2 , but an absurd energy, lower than E%. Concavity is a 
mandatory condition for the universality of T . 

One way to simultaneously implement pairing corrections and use a mini- 
mal paraboloid term consists in adding to the physical H terms such as, p x 
[cos 2 (tt N/2) + cos 2 (7rZ/2)-l]+«N 2 + 2uNZ + wZ 2 . (Notice how the cosines 
take advantage of the continuous nature of (N) and (Z) at finite temperatures.) 
Concavity is achieved if u,v,w are tuned to compensate for the worst tuned 
SD minima listed above, namely { — 1500, —930, —980, —120}. A rough solution 
consists in taking p ~ 1100 and a "circular" ansatz, P x (N 2 + Z 2 ) , with 
P ~ 800, but less brutal solutions are obviously worthy of consideration, in 
order to minimize this tampering with the Hamiltonian. In any case, we stress 
that the operators, N and Z, commute with H, and, therefore, that the present 
counterterms do not perturb nuclear dynamics. Once again we emphasize that 
this approach unifies the treatment of odd and even nuclei. The fact that the 
nuclear density functional needs to originate from a Hamiltonian, completed by 
such counterterms, illustrates how important concavity is for making predictions 
for nuclei far from stability, i.e., for the so-called exotic nuclei. 

5 Summary, discussion and conclusion 

We have demonstrated how a list of ground-state energies for a sequence of 
isotopes can be turned into a concave pattern. This involves simple manipula- 
tions; for instance, an explicit term, accounting for pairing in even nuclei, can 
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be subtracted from the bindings. This unifies the treatment of odd and even 
nuclei, a notoriously difficult problem. If needed, a small quadratic correction 
can also be added to guarantee concavity at all points along the sequence. 

Similar arguments leading to concavity clearly hold for isotones as well, 
and, furthermore, for any other sequence of neighboring nuclei in any direction 
across the nuclear table. Once this empirical tuning has been implemented, 
linear (or more general) extra- and interpolations of the concave pattern can 
provide surprisingly accurate and robust estimates of, or bounds for, binding 
energies. These tuning terms, which are added to induce concavity, are, of 
course, subtracted in fine. 

This work then defined a more ambitious extra- and interpolation scheme, 
involving thermodynamical functions from a grand canonical ensemble, because 
such functions may have rigorous concavity properties. Theorems are, indeed, 
available to prove such properties. For instance, the free energy is a concave 
function of the average particle number and is also a decreasing function of the 
temperature. We also found strong numerical evidences concerning the average 
energy. This average energy at nonzero temperature turns out, in general, to be 
a concave function of the average particle number, except for "minor accidents," 
which depend upon the temperature. 

For every given, finite temperature, we found that the average energy and 
the free energy, as functions of the average particle number, give upper and 
lower bounds, respectively, for the concave envelope of the ground-state energies. 
When the temperature vanishes, both bounds converge to the exact results. At 
this vanishing temperature, however, the analyticity of such thermodynamical 
functions is lost, because their limit is only piecewise continuous. It is, therefore, 
necessary to retain a minimum temperature if one wants to obtain practical 
extrapolations for the prediction of exotic nuclei. A minimum amount of particle 
number fluctuation is necessary to validate the conversion of particle number, 
an initially discrete quantity, into a continuous variable. 

We, therefore, implemented numerical estimates of several thermodynamical 
functions at moderate temperatures, a few hundred keV at most. This yields 
a first result, namely a "band", enclosing ground-state energies between the 
average energy and the free energy. The width of the band defines an error bar 
which can be trusted when extrapolations are made. 

A difficulty arises, however, because of an insufficient number of excited 
states for nuclei at both ends of any sequence of isotopes. These excited-state 
energies are essential for calculating a meaningful value of the average energy 
for any value of A. More often, only the ground-state energy is known for such 
neutron very rich or neutron very poor nuclei. The calculation of the average 
and free energies is, thus, possible only in an interval smaller than the inter- 
val of masses where ground states are known. Two tactics are then available, 
namely i) an extrapolation only of the sequence of ground-state energies and 
ii) an extrapolation of the thermodynamical functions, starting from a smaller 
interval. We found that the first tactic, namely at zero temperature, is reliable, 
especially when completed by the error bar derived from the finite temperature 
considerations. This work also tested the second tactic, with some success but 
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with the limitation due to the lack of known excited states. The lack leads 
to edge effects in the calculation of thermodynamical functions. Therefore, if 
linear, or polynomial, or more general analytical fits of such thermodynamical 
functions are attempted, it is better to restrict them from a slightly smaller 
mass interval. Then, extrapolations should be reliable within one and maybe 
two mass units towards drip lines. 

We can make the strong conclusion that the combination of concavity and 
extrapolations of thermodynamical functions gives a systematic set of upper 
and lower bounds for the prediction of ground-state energies. 

Our last, and perhaps most important result, is the connection between 
concavity and the universality of the density functional. In a two-dimensional 
analysis, we showed that a paraboloid term with respect to the proton and 
neutron number operators must be added to the Hamiltonian, so as to guarantee 
obtaining consistent energy minima everywhere in density space. This term can 
be made minimal if counterterms for pairing are also added. 
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